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We report a measurement of the branching fractions for B — > (-7r)£~^ decays based on 341.1 
fb _1 of data collected at the T(4S) resonance with the BABAR detector at the PEP-II e + e~ storage 
rings. Events are tagged by fully reconstructing one of the B mesons in a hadronic decay mode. 
We obtain B(B- -f D°ij> e ) = (2.33 ± 0.09 stat . ± 0.09 syst .)%, 6(B~ ->■ D*H~vi) = (5.83 ± 
0.15 stot . ±0.30 syst .)%, B{B° -> D+e-ut) = (2.21 ± 0.1W ± 0.12 syst .)%, B(B° -> D*+£-p e ) = 
(5.49 ± 0.16 stot . ±0.25 syst .)%, B(B~ D + i:-i-v t ) = (0.42 ± 0.06 siat . ± 0.03 syst .)%, B(B~ 
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D*+ir-e-Pt) = (0.59±0.05 stat . ±0.04 syBt .)%, B(B° -> D°Ti + rv e ) = (0.43 ± 0.08 stat . ± 0.03 sys 
and B(B° -> D*°Tv + ru e ) = (0.48 ± 0.08 stat . ± 0M syst .)%. 

PACS numbers: 13.20He,12.38.Qk,14.40Nd 



The determination of the individual exclusive branch- 
ing fractions of B — ► X c £~Vi decays [l[ is important for 
the study of the scmilcptonic decays of the B meson. 
Improvement in the knowledge of these branching frac- 
tions is also important to reduce the systematic uncer- 
tainty in the measurements of the Cabibbo-Kobayashi- 
Maskawa matrix elements \V c b\ and \V u b\- For ex- 
ample, one of the leading sources of systematic uncer- 
tainty in the extraction of \V c b\ from the exclusive de- 
cay B — > D*l~vi is the limited knowledge of the back- 
ground due to B — > D*ir£~Di. Improved measurements 
of B — > X c t~i>i decays will also benefit the accuracy of 
the extraction of \V u b\, as analyses are extending into 
kinematic regions in which these decays represent a siz- 
able background. 

Based on current measurements 0, 0, H, H, 0] the rate 
of inclusive semileptonic B decays exceeds the sum of the 
measured exclusive decay rates [8]. While B — » D£~Dg 
and D*l~i>£ decays account for about 70% of this total, 
the contribution of other states, including resonant and 
non-resonant D^n£~v t (denoted by D**£~vt), is not 
yet well measured and may help to explain the inclusive- 
exclusive discrepancy. 

In this letter, we present measurements of the branch- 
ing fractions for B — > D^(n)£~Di decays, separately for 
charged and neutral B mesons. 

The analysis is based on data collected with the BABAR 
detector Q at the PEP-II asymmetric-energy e + e~ stor- 
age rings. The data consist of a total of 341.1 fb _1 
recorded at the T(4S) resonance, corresponding to 378 
million BB pairs. An additional 36 fb _1 off-peak data 
sample, taken at a center-of-mass (CM) energy 40 MeV 
below the T(AS) resonance, is used to study background 
from e + e~ — > // (/ = u,d,s,c,T) events (continuum 
production). A detailed GEANT4-based Monte Carlo 
(MC) simulation of BB and continuum events is 
used to study the detector response, its acceptance, and 
to test the analysis techniques. The simulation models 
B^D^l-Di decays using calculations based on Heavy 
Quark Effective Theory [ll|, B — > £)**(-> D^ir)e-p e 
decays using the ISGW2 model [l2j , and B — > D^TTl~v t 



decays using the Goity-Roberts model [lc 

We select semileptonic B decays in events containing 
a fully reconstructed B meson (B tag ), which allows us to 
constrain the kinematics, reduce the combinatorial back- 
ground, and determine the charge and flavor of the signal 
B meson. 

We first reconstruct the semileptonic B decay, selecting 
a lepton with momentum p| in the center-of-mass frame 
higher than 0.6 GeV/c. Electrons from photon conver- 
sions and 7T° Dalitz decays are removed by searching for 



pairs of oppositely charged tracks that form a vertex with 
an invariant mass compatible with a photon conversion or 
a 7T° Dalitz decay. Candidate D° mesons, having the cor- 
rect charge-flavor correlation with the lepton, are recon- 
structed in the K~ir + , if~7r + 7r°, K ~Tr + ir + ir~ , i^c7r + 7r _ , 



i^TT+TT-TT , K° s lt a , K+K- 



and K S K S channels, 



TT+TT+TT , K° S TT+, 



7T 7T 
, 7T+7T 

and D + mesons in the K~-k + tt + , K 
K s ir+ir , K+K-TT+, Kp[ + , and ^tt+tt+tt" channels. 
In events with multiple B — » D£~9i candidates, the can- 
didate with the best D-£ vertex fit is selected. Candidate 
D* mesons are reconstructed by combining a D candi- 
date with a pion or a photon in the D* + — ► D°ir + , 
D*+ -► L>+7T°, D*° -> D°7r°,_and D*° -> D° 1 chan- 
nels. In events with multiple B — > D*£~vi candidates, 
we choose the candidate with the smallest % 2 based on 
the deviations from the nominal values of the D invariant 
mass and the invariant mass difference between the D* 
and the D, using the measured resolution. 

We reconstruct B tag decays of the type B — ► DY, 
where Y represents a collection of hadrons with a total 
charge of ±1, composed of rii7r ± + n 2 K ± +n 3 K s + 71471"°, 
where ni + n 2 < 5, n 3 < 2, and n± < 2. Using D°(D + ) 
and D*°(D* + ) as seeds for B~(B°) decays, we recon- 
struct about 1000 different decay chains. 

The kinematic consistency of a B tag candidate with 
a B meson decay is evaluated using two variables: the 
beam-energy substituted mass tues = V / s/4 — bspj 
and the energy difference AE — E* B — \fs/2. Here 
•^/i refers to the total CM energy, and p* B and E* B de- 
note the momentum and energy of the B tag candidate 
in the CM frame. For correctly identified B tag decays, 
the niES distribution peaks at the B meson mass, while 
AE is consistent with zero. We select a B tag candidate 
in the signal region defined as 5.27 GeV/c 2 < tiies < 
5.29 GeV/c 2 , excluding B tag candidates with daughter 
particles in common with the charm meson or the lepton 
from the semileptonic B decay. In the case of multiple 
B tag candidates in an event, we select the one with the 
smallest \AE\ value. The B tag and the D^*'£ candidates 
are required to have the correct charge-flavor correlation. 
Mixing effects in the B° sample are accounted for as de- 
scribed in [14fl . Cross-feed effects, i.e., B^ ag (B^ ag ) can- 



didates erroneously reconstructed as a neutral (charged) 
B, are subtracted using estimates from the simulation. 

For B — > D^X£~Di decays, D(D*) candidates are 
selected within 2a (1.5-2.5(7, depending on the D* decay 
mode) of the D mass (D* — D mass difference), with a 
typically around 8 (1-7) MeV/c 2 . We also require the 
cosine of the angle between the directions of the 
candidate and the lepton in the CM frame to be less 
than zero, to reduce background from non-B semileptonic 
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decays. 

We reconstruct B~ -> D^+ir~£-i/ e and B° -> 
n + £~ Di decays starting from the corresponding 
B — ► D^*'X£~p£ samples and selecting events with only 
one additional reconstructed charged track that has not 
been used for the reconstruction of the B tag , the sig- 
nal £>(*), or the lepton. For the B Q -> D°n + £~vi 



and the B° -> D* 



Kg decays, we additionally re- 



quire the invariant mass difference M(Dn) — M(D) to 
be greater than 0.18 GeV/c 2 to veto B° -> D* + £-p e 
events. To reduce the combinatorial background in the 
B° — > D*°ir + £~ise mode, we also require the total extra 
energy in the event, obtained by summing the energy of 
all the showers in the electromagnetic calorimeter that 
have not been assigned to the B tag or the D^*>£" candi- 
dates, to be less than 1 GeV. 

The exclusive semileptonic B decays are identified 
by the missing mass squared in the event, m 2 niss — 
(p(T(45)) - p(B tag ) - p(£)W(7r)) - pii)) 2 , defined in 
terms of the particle four-momenta in the CM frame 
of the reconstructed final states. For correctly recon- 
structed signal events, the only missing particle is the 
neutrino, and rn 2 niss peaks at zero. Other B semileptonic 
decays, where one particle is not reconstructed (feed- 
down) or is erroneously added (feed-up) to the charm 
candidate, exhibit higher or lower values in m 2 niss . To 
obtain the B semileptonic signal yields, we perform a 
one-dimensional extended binned maximum likelihood 
fit [15(] to the rr? misl . distributions. The fitted data sam- 



ples are assumed to contain four different types of events: 
B — > D^*\-k)£~Di signal events, feed-down or feed-up 
from other B semileptonic decays, combinatoric BB and 
continuum background, and hadronic B decays (mainly 
due to hadrons misidentified as leptons). For the fit to the 
m miss distributions of the B — > D^n£^De channel, we 
also include a component corresponding to other misre- 
constructed B — > D** (D*ir)£~ Di decays. We use the MC 
predictions for the different B semileptonic decay rr? mil , s 
distributions to obtain the Probability Density Functions 
(PDFs). The combinatoric BB and continuum back- 
ground shape is also estimated by the MC simulation, 
and we use the off-peak data to provide the continuum 
background normalization. The shape of the continuum 
background distribution predicted by the MC simulation 
is consistent with that obtained from the off-peak data. 

The m 2 niss distributions are compared with the results 
of the fits in Fig. 1 for each of the B -> (ir)£- v e 
channels. The fitted signal yields and the signal efficien- 
cies, accounting for the B ta g reconstruction, are listed in 
Table H 

To reduce the systematic uncertainty, the exclusive 
B(B — ► (n)£~ Pi) branching fractions relative to the 
inclusive semileptonic branching fraction are measured. 
A sample of B — > X£~p£ events is selected by identifying 
a charged lepton with CM momentum greater than 0.6 
GeV/c and the correct charge- flavor correlation with the 
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FIG. 

B~ D { '£ 
B° -> D* + £~vi 

R) B° D°TT 



(Color online) Fit to the m 

,0e--.. b ) B - ^ D*°£~ve, 
e) B~ — » D + ir~£~v t 
£-u e , and h) B° 



2 

mis 

c) 



distribution for a) 
B° -)• D+Tue, d) 
f) B~ ~> D*+-K-£-v t , 
D*°iv + rvf. the data 



(points with error bars) are compared to the results of the 
overall fit (sum of the solid histograms). The PDFs for the 
different fit components are stacked and shown in different 
colors. 



B tag candidate. 



In the case of multiple B tag candidates 
in an event, we select the one reconstructed in the decay 
channel with the highest purity, defined as the fraction of 
signal events in the tues signal region. Background com- 
ponents peaking in the tubs signal region include cascade 
B meson decays [i.e., the lepton does not come directly 
from the B) and hadronic decays, and are subtracted 
by using the corresponding MC distributions. The total 
yield for the inclusive B — » X£~vi decays is obtained 
from a maximum likelihood fit to the tues distribution 
of the B tag candidates using an ARGUS function [l6| 
for the description of the combinatorial BB and contin- 
uum background, and a Crystal Ball function [l7| for the 
signal. Additional Crystal Ball and ARGUS functions 
are used to model a broad-peaking component, included 
in the signal definition, due to real B — > X£~Di decays 
for which, in the B tag reconstruction, neutral particles 
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have not been identified or have been interchanged with 
the semileptonic decays. Fig. [2] shows the tues distri- 
bution of the B tag candidates in the B~ — * X£~T>i and 
W -> Xl~ v i sample. The fit yields 159896 ± 1361 events 
for the B~ — > X£~vi sample and 96771 ± 968 events for 
the B° — > Xl~i>e sample. 




B(B — » Xl~vi) used for the normalization (1.9%). We 
also include a systematic uncertainty due to differences 
in the efficiency of the B tag selection in the exclusive 
selection of B — > D^*\t:)£~V£ decays and the inclusive 
B — ► Xl~vi reconstruction (0.9-5.6%), and the extrac- 
tion of the B -► D^{-K)l-v t (0.4-1.8%) and B -> Xl~v t 
(0.5-0.9%) signal yields. The complete set of systematic 
uncertainties is given in Ref . [18( . 

We measure the following branching fractions: 



FIG. 2: (Color online) ities distributions of the a) B~ — » 
Xi~vi, and b) _B° — » XtTDt samples. The data (points with 
error bars) are compared to the result of the fit (solid line). 
The dashed lines show the broad-peaking component and the 
sum of the combinatorial and continuum background. 



The relative branching fractions B(B — > 
D(*'(jr)l~Dt)/B(B — > Xt~Di) are obtained by correcting 
the signal yields for the reconstruction efficiencies 
(estimated from BB MC events) and normalizing to the 
inclusive B — > X£~V£ signal yield, following the relation 
B(B -> DW{-k)1-Vi)/B{B -> Xl-vt) = Here, 

N*i„ is the number of B -> L> ( *) (tt)£" ' signal events for 



B(B 


- D°r 


~Vl) 


= (2.33±0.09 stat . 


±0.09 syst .) 


% 


B(Br 


-> z>*°r 


~vi) 


= (5.83±0.15 stat . 


±0.30 syst .) 


% 


B(B° -> £> + r 


~Vl) 


= (2.21 ± 0.1W 


± 0.12 SJ/S t. 


% 


B{B° 


-> zr+r 


~Vl) 


= (5.49±0.16 stai . 


± 0.258^84. 


% 


B{B- - 




~vi) 


= (0.42±0.06 siat . 


±0.03 SJ/st .) 


% 


B{B~ -> 




~Vl) 


= (0.59±0.05 stat . 


±0.04 SJ/St .) 


% 


B(B°- 




~vt) 


= (0.43±0.08 stat . 


±0.03 syst .) 


% 


B(B°^ 


D* n+t 




= (0.48±0.08 s t at . 


±0.04 syst .) 


% 



The accuracy of the branching fraction measurements 
for the B — > D^l~vi decays is comparable to that 
of the current world average ||. We compute the to- 
tal branching fractions of the B — > D^itfTVi decays 
assuming isospin symmetry, BiB — ► D^*'tt°£~U£) — 



the various modes, reported in Table [J together with the \ B{ ? ~* ^ estimate the branching frac 



corresponding reconstruction efficiencies e s 



is the tions of D^*'n° final states, obtaining: 



B XI i/£ signal yield, and e s i is the corresponding 
reconstruction efficiency including the B tag reconstruc- 
tion, equal to 0.36% and 0.23% for the B~ -> Xl~v t 
and B° — > Xl~vi decays, respectively. The absolute 
branching fractions B(B — ► D^*\tt)£~Di) are then 
determined using the semileptonic branching fraction 
B(B -► A£-p £ ) = (10.78 ± 0.18)% and the ratio of the 
B° and the 5+ lifetimes t b + /t b o = 1.071 ± 0.009 [§]. 

Numerous sources of systematic uncertainties have 
been investigated. The uncertainties due to the detec- 
tor simulation are established by varying, within bounds 
given by data control samples, the tracking efficiency of 
all charged tracks (resulting in 1.2-2.7% relative system- 
atic uncertainty among the different decay modes), the 
calorimeter efficiency (0.5-1.8%), the lepton identification 
efficiency (0.4-3%), and the reconstruction efficiency for 
low momentum charged (1.2%) and neutral pions (1.3%). 
We evaluate the systematic uncertainties associated with 
the MC simulation of various signal and background pro- 
cesses: photon conversion and ir° Dalitz decay (0.04- 
0.4%), B cascade decay contamination (0.6-1%), and fla- 
vor cross-feed (0.2-0.3%). We vary the B — > D£~vi and 
B — > D*£~p£ form factors within their measured uncer- 
tainties [HI (0.4-0.8%) and we include the uncertainty on 
the branching fractions of the reconstructed D and D* 
modes (2.3-4.4%), and on the absolute branching fraction 



B(B~ 
B(B° 



■ D^irl've) = (1.52±0.12 fltat . ±0.10 fl „ at .)% 
>DV>wtTvt) = (1.37±0.17 stat . ±0.10 s 



syst.)%, 



where we assume the systematic uncertainties on the 
B — > Dirt~i>i and B — > D*n£~Pi modes to be completely 
correlated. These results are consistent with, but have 
smaller uncertainties than, recent results from the Belle 
collaboration 

By comparing the sum of the measured branching frac- 
tions for B — > (tt)£~ vi with the inclusive B — > 
X C £~V>£ branching fraction [§], a (11 ± 4)% discrepancy 
is observed, which is most likely due to B — > D^nitfTui 
decays with n > 1. 
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TABLE I: Signal yields and reconstruction efficiencies for the B — > D^(n)£ vg decays. 



Decay Mode 


N Sl g 


e sl9 (xHT 4 ) 


Decay Mode 


iV 3 ; 9 


e al9 (xHT 4 ) 


B~ -> D u n? e 


f635 ± 6f 


f .7f ± 0.02 


B~ -> D + TY~rui 


f 74 ± 25 


f .02 ± 0.03 


B~ -> D M Tv t 


3050 ± 73 


f .27 ± O.Of 


B~ -» D* + ty-£-ui 


306 ± 27 


f .26 ± 0.03 


B° -> D+£~vi 


852 ± 40 


0.94 ± 0.02 


B u D u ir+£-v e 


f 07 ± 20 


0.60 ± 0.03 


B° -> D*+£~v t 


2045 ± 55 


0.9f ± O.Of 




f30 ± 20 


0.66 ± 0.02 
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TABLE II: Systematic uncertainties (relative errors in %) in the measurement of B(B -> 't~ v t )jB(B -> Xi~v t ). 





Systematic uncertainty on B{B 


-> D [ -* ) i-v l )IB(B -> Xl~v t ) 




B~ -> D rvt 


B~ -> D*°rv e 


B° -> D+e-ue 


B° -> D* + rv t 


Tracking efficiency 


1.4 


1.2 


1.4 


1.5 


Neutral reconstruction 


0.7 


1.9 


0.5 


1.1 


lepton ID 


0.5 


0.4 


0.5 


0.6 


Soft particle efficiency 




1.3 




1.2 


Monte Carlo corrections 




Conversion and Dalitz decay background 


0.04 


0.07 


0.06 


0.05 


Cascade B — > X — > i~ decay background 


0.6 


0.6 


1.0 


1.0 


B° — B~ cross- feed 


0.2 


0.3 


0.2 


0.3 


Form factors 


0.4 


0.8 


0.4 


0.8 


D branching fractions 


2.3 


2.1 


4.1 


2.6 


D* branching fractions 




2.8 




0.8 


B — > Xl~vi branching fraction 


1.9 


1.9 


1.9 


1.9 


B tag selection 


0.9 


1.7 


1.8 


1.3 


Fit technique 




B -> Xl'vt yield 


0.5 


0.5 


0.9 


0.9 


B -> D M rv e yield 


0.6 


0.4 


1.2 


0.4 


Total systematic error 


3.7 


5.2 


5.4 


4.5 



TABLE III: Systematic uncertainties (relative errors in %) in the measurement of B(B — » D^nl vt)/B{B — * XI vt). 





Systematic uncertainty on B{B - 


-» D(*Ke-& e )/B(B - 






b~ D + TT~e~p e 


B- -^D* + TT-ri> e 


B» -> D°TT + rv e 


b° -» L>* u 7r+r^ 


Tracking efficiency 


1.8 


2.7 


1.5 


1.7 


Neutral reconstruction 


1.7 


1.8 


1.1 


1.8 


lepton ID 


2.3 


3.0 


2.6 


1.8 


Soft particle efficiency 




1.2 




1.3 


Monte Carlo corrections 




Conversion and Dalitz decay background 


0.15 


0.4 


0.05 


0.2 


Cascade B — > X — > £~ decay background 


0.6 


0.6 


1.0 


1.0 


B° — B~ cross- feed 


0.2 


0.3 


0.2 


0.3 


Form factors 


0.4 


0.8 


0.4 


0.8 


D branching fractions 


4.2 


2.8 


2.5 


2.9 


D* branching fractions 




0.9 




3.3 


B — » Xl~vi branching fraction 


1.9 


1.9 


1.9 


1.9 


Btag selection 


5.0 


4.3 


4.0 


5.6 


Fit technique 




B -> Xl~v t yield 


0.5 


0.5 


0.9 


0.9 


B -> D { *Krv t yield 


1.2 


0.9 


1.8 


1.5 


Total systematic error 


7.7 


7.3 


6.4 


8.4 



